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ABSTRACT 



We present Spitzer Space Telescope time series photometry of the exoplanet system 
HD 189733 spanning two times of secondary eclipse, when the planet passes out of view 
behind the parent star. We estimate the relative eclipse depth in 5 distinct bands and 
find the planet-to-star flux ratio to be 0.256±0.014% (3.6 //m), 0.214±0.020% (4.5 fim), 
0.310 ± 0.034% (5.8 //m), 0.391 ± 0.022% (8.0 //m), and 0.598 ± 0.038% (24 //m). For 
consistency, we re-analyze a previously published time series to deduce a contrast ratio 
in an additional band, 0.519 it 0.020% (16 //m). Our data are strongly inconsistent with 
a Planck spectrum, and we clearly detect emission near 4 //m as predicted by published 
theoretical models in which this feature arises from a corresponding opacity window. 
Unlike recent results for the exoplanet HD 209458b, we find that the emergent spec- 
trum from HD 189733b is best matched by models that do not include an atmospheric 
temperature inversion. Taken together, these two studies provide initial observational 
support for the idea that hot Jupiter atmospheres diverge into two classes, in which a 
thermal inversion layer is present for the more strongly irradiated objects. 

Subject headings: binaries: eclipsing — infrared: stars — planetary systems — stars: 
individual (HD189733) — techniques: photometric 



1. Introduction 



Of t he many intriguing avenu es for investigation afforded by transiting exoplanets (for a re- 



view, see 



Charbonneau et al.ll2007l ). perhaps the most exciting is the opportunity to study directly 
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the planetary atmosphere without the need to spatially resolve the light from the planet from 
that of its host star. At visible wavelengths, where the emitted light from the planet is negligi- 
ble, the technique of transmission spectroscopy has been employed from space-based platforms to 



study atoms, molecules, and clouds in th e atmospheres of H P 209458b (ICharborineau et al 



Vidal-Madiar et al.1 l2003l : iBarmanI 12003) and HD 189733b (jTinetti et al 



2007 



Pont et alj 



2002 : 



200S : 



dMoutou et alJ 


2001, 


2003: 


Winn et al. 


2004 


Deming et al. 


2005a: 


Narita et al. 


2005: 


Arribas et al. 



Redfield et al 



(2003) of 



200d : iBozorgnia et al.l bood ) until the recent breakthrough detection by | 
atomic sodium in the atmosphere of HD 189733b. Searches for reflected starlight have deliv- 
ered incr easingly stringent constra i nts on the wavelength-dependent geometric alb ed o and phase 



2006 



function (ICharbonneau et al.l Il999l : ICollier Cameron et al. 



2002 



Leigh et al. 



2003 



iRowe et al 



20081 ). At infrared wavelengths, the planet-to-star contrast ratio improves dramatically, per- 
mitting a direct study of the emitted planetary radiation through the modulation of the total 
light from the star and planet spanning times of secondary eclipse, when the planet is occulted 
by the star. Ground-based sear ch es for this eff e ct have been frustrated to date by telluric vari- 



ability (IRichardson et al 



2003aH bl: ISnellenI l200,4 IPeming et al.l l2007al : iKnutson et al.l l2007bl ^. al 



though ISnellen &: Covind (120071 ) have recently claimed a tentative detection for 0GLE-TR-113b. 
Exem plifying the power o f a small-aperture telescope in an ideal location, the Spitzer Space Tele- 
scope ([Werner et al.l |200J) has revol utionized this field by yielding a flurry of det e ctions of the 



emitted radiation from 6 exoplanets (ICharbcinnea.u et al 



Demorv et al.l 120071 : [Harrington et al 



200i 



determ inations of the emitted spectra (jGrillmair et al 



2007 



2005; 



Knutson et al 



2007 



2007a 



Deming et al.ll2005bl . l2006l . l2007bl : 



2008 ) , including recent 



Richardson et al.l 120071 : ISwain et al 



2008bl ) of HD 189733b and HD 209458b. Numerous ongoing programs promise additional exciting 
results to come. Although Spitzer will exhaust its cryogen in the spring of 2009, two photometric 
channels will continue to operate at full strength and w ill provide a singular window into the at- 
mospheres of these distant worlds (iDeming et al.ll2007d ). if the warm phase of the Spitzer Mission 
is approved. 

Each of the techniques described above is facilitated for large planets that orbit nearby, bright 
stars at small orbital separations. The discovery of the exoplanet HD 189733b (jBouchv et al.ll2005l ) 
provided a golden opportunity for practitioners of such techniques for two reasons. First, it is at a 
distance of only 19.3 pc, and hence at infrared wavelengths it is the brightest star (K = 5.5) known 
to host a transiting exoplanet. Second, the planet-to-star surface area ratio is th e second largest of 



the 2 8 transiting exoplanets described in the literature, surpassed by only TrES-3 (lO 'Donovan et al 



20071 ). It is primarily for these reasons of accessibility that we selected HD 189733b for intensive 
study with Spitzer. For the purposes of this paper, we adopt the values for orbital inclination 
i, the planet radi u s Rp, the stellar radius R^,, and the time of transit center as deter mined 
bv iKnutson et al.l (|2007al ). and the orbital period P as determined bv IWinn et al.l (j2007al ). We 
note that the targ et star has an M-dwarf companion at a projected physical separation of 216 AU 
(jBakos et al.ll2006l ). This companion appears in the field of view of the data we present below, and 
we comment upon its impact upon our photometric results when appropriate. 
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The recent Spitzer detections of the flux emitted by several exoplanets have spurred numerous 
theoretical efforts t o model their atmospheres and predict their emitted spectra (see for ex ample 



)sp 

i, 



Barman et alJl2005l : Burrows et alJl2005l . 12009 : iFortnev et al.ll2005l . l2006l : ISeaeer et al.ll2005l l. The 



constraints on the theoretical predictions become much more stringent with detections at numerous 
wavelengths, from which a broadband spectrum may be inferred. In this paper, we describe our 
estimate of planet-to-star contrast ratio in bands centered at 3.6, 4.5, 5.8, 8.0, and 24 um. We then 
combine these results with our re-analysis of the earlier measurement at 16 fim (iDeming et al.ll2006l ) 
to consider what constraints may be placed on the temperature and chemistry of the planetary 
atmosphere, and the efficiency with which it redistributes the incident stellar radiation from the 
da y side to the n i ght s ide. These observations are particularly timely, given the recent study 
bv iKnutson et al.l (j2008l ). which included data presented in Deming et al. (2005b), of the broad 
band spectrum of HD 209458b in 5 of the above bandpasses. In the assembled spectrum, features 
that had been predict ed in absorption were observed in ern i ssion , a cl ear sign of an a tmosp heric 
temperature inversion ([Burrows et al.ll2007al ). iFortnev et al.l (|2008l ) and lBurrows et al.l (|2008l ) have 
recently proposed a division between the predicted atmospheric temperature-pressure profiles for 
hot Jupiters dependent upon the level of stellar insolation, in which only the m ost strongly irradiate d 
planets show a temperature inversion. According to the criterion proposed by lFortney et al.l (|2008l ). 
HD 209458b is a member of the pM class and hence is expected to show a temperature inversion, 
whereas HD 189733b is a member of the pL class and should not. Intriguingly, HD 209 458b and 
HD 189733b also differ in that the radi us of the former (1.32 .Rjnp; iKnutson et al.ll2007d ) is much 
larger than that of the latter (1.14i?jup; Knutson et al.ll2007al ) and is at odds with the predictions 
of models that do not include an additional energy source interior to the planet. Whether the large 
radius and the existence of a temperature inversion share a common physical origin is unclear, 
but surely worthy of further investigation. Our study of the broad band spectrum of HD 189733b 
provides an ideal opportunity to conduct a comparative study with HD 209458b. 



2. Observations and Analysis 
2.1. IRAC Observations 



We observed HD 189733b over a period of 4.9 hours on UT 2005 November 24, spanning a 
single secondary eclipse, using the Infrared Array Camera (IRAC; iFazio et al.ll2004l ) on the Spitzer 
Space Telescope. We observed in subarray mode with an exposure time of 0.1 s and cycled between 
the four IRAC channels in order to obtain estimates of the depth of the eclipse at 3.6, 4.5, 5.8, 
and 8.0 /im simultaneously. We obtained a total of 21,760 32 x 32 pixel images in each channe0. 



We use images processed using version S13.0 of the standard Spitzer pipeline, to avoid the additional noise 
introduced by a dark drift correction that was applied to all subarray images beginning with version S14.0 of the 
pipeline, released in May 2006. This dark drift correction is poorly constrained for subarray images dominated by a 
single bright star, and as a result introduces noise at a level higher than the effect it is meant to correct. See IRAC 
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Images are taken in sets of 64, and four sets of images (4 x 64 images total) are obtained in each 
channel before re-pointing the telescope to position the star on the subarray for the next channel. 
The position of the star on the subarray is still varying through the first set of 64 images after 
each re-pointing, and we chose to discard this initial set of images, leaving a total of 16,320 usable 
images in each of the four channels. The total size of this pointing drift is 0.3 pixels in the first 
set of 64 images, and 0.1 pixels or less in the following three sets of images. For each image, we 
calculate the JD value for the time at mid-exposure, and apply a correction to convert these JD 
values to the appropriate HJD taking into account the position of Spitzer in the Solar system at 
each point during the observations. 

Because the two shortest wavelength IRAC channels (3.6 and 4.5 fim) use InSb detectors and 
the two longer wavelength channels (5.8 and 8.0 fim) use Si:As detectors, there are fundamental 
differences between the properties of the data taken with these two types of detectors. We describe 
our analysis for each type of detector separ ately below. The specific methods we employ our similar 



to those described in iKnutson et al.l (j2008l ). 



2.1.1. 3.6 and 4.5 Observations (InSb Detector) 



Because HD 189733 is a bright star and the background at these shorter wavelengths is minimal, 
we calculate the flux from the star in each image using aperture photometry with a circular aperture 
with a radius of five pixels. We determine the position of the star in each image as the position- 
weighted sum of the flux in a 7 x 7 pixel box centered on the approximate position of the star. We 
estimate the background in each image by selecting a subset of pixels from the corners of the image 
where the point spread function of the star is faintest and avoiding the region around the M dwarf 
companion, making a histogram of the flux values in these pixels, and fitting a Gaussian function 
to the center of this distribution. 

Fluxes measured at these two wavelengths show a strong correlation with the changing position 
of the star on the array, at a level comparable to th e depth of the sec o ndary eclipse. This effect 
is due to a well-docuniented intra-pixel seiisitivit y (jReach et al.l l2005l : ICharbonneau et al.l l2005l : 
Morales-Calderon et al.ll2006l : iKnutson et al.ll2008l ). and can be removed by fitting the data with a 
quadratic function of the sub-pixel position of the center of light of the stellar image, according to: 

= / [ci + C2 (x - 13.5) + C3 (x - 13.5)2 + C4 (y - 14.5) + C5 (y - 14.5)^] (1) 

where / is the flux from the star prior to this distortion, is the measured flux, x and y denote the 



pipeline history available online at littp://ssc. spitzer. caltech.edu/irac/ for more information. This is important for 
only the 3.6 /im channel, where the eclipse depth is partially degenerate with the decorrelation as a function of x and 
y position on the array. If we use the new version of the pipeline with the poorly constrained dark drift correction, 
we find the measured x and y positions of the star on the array are shifted by up to 0.001 pixels and the depth of 
the eclipse in this channel increases to 0.0041. 
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location in pixels of the center of the star on the subarray, and ci 
in the fit. 



C5 are the five free parameters 



We fit the correction for the intra-pixel sensitivity of the array and the eclipse c urve simul- 
taneously to the d ata using a Markov Chain Monte Carlo method (see for example iFordI l2005l : 
Winn et alJl2007bl ) with 10^ steps. We set the uncertainty on individual points equal to the stan- 
dard deviation of the out-of-eclipse data after correction for the intra-pixel variations, and remove 
outliers of 5a or more as calculated using the residuals from the best-fit light curve. We allow both 
the depth and timing of the secondary eclipse to vary independently for the eclipse s at each of the 
two ob served wavelengths, and fix the v alues of i, Rp , R^,, a nd to those given in lKnutson et al 



(|2007al ). the value f or P to that found by 



Winn et al 



(|2007al ). We calculate our eclipse curve using 



the equations from lMandel &: Agoll (j2002l ) for the case with no limb-darkening. After running the 
chain, we search for the point in the chain where the value first falls below the median of all 
the values in the chain (i.e. where the code had first found an excellent fit), and discard all the 
steps up to that point. We take the median of the remaining distribution as our best-fit parameter, 
with errors calculated as the symmetric range about the median containing 68% of the points in 
the distribution. The distribution of values was very close to symmetric in all cases, and there 
did not appear to be any strong correlations between variables. Figure [1] shows the binned data 
with the best fit to the detector effects overplotted, and Figure [2] shows the binned data once these 
trends are removed, with best-fit eclipse curves overplotted. The best-fit eclipse depths and times 
are given in Tabled) 



2.1.2. 5.8 and 8.0 jjLm Observations (Si:As Detector) 



At longer wavelengths the fiux from the star is smaller and the zodiacal background is larger; 
as a result we chose to use a smaller circular aperture with a radius of 3.5 pixels in order to minimize 
the noise contribution from this increased background. As before, we calculate the position of the 
star individually in each image as the position-weighted sum of the fluxes in a 7 x 7 pixel box, and 
estimate the background using a Gaussian fit to a histogram of the pixels in the corners of the 
array. Fluxes in the first 10 images and the 58th image in each set of 64 are consistently below the 
median value for the set by as much as 10%, with the lowest values at the beginning of each set, 
so we chose to exclude these 11 images from each set of 64 in our analysis. 



There is no known intra-pixel sensitivity at t hese wavelengths (jCharbonneau et al.ll2005l ). but 
there is another well-documented detector effect (jKnutson et al.ll2007al ). which causes the effective 
gain (and thus the measured fiux) in individual pixels to increase over time. T his effect has 



been referred to as the "detector ramp", and has also been observed ([Peming et al.l 120061 ) in the 
IRS 16 (j,m. array, which is made from the same material. The size of this effect depends on the 
illumination level of the individual pixel, such that the time-dependent gain of different pixels form 
a well-defined family of functions. The most strongly-illuminated pixels (>250 MJy Sr^^ in the 
8 /im channel) converge asymptotically to a constant value within the first hour of observations. 
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The time scale for this convergence increases with decreasing illumination, and the limiting case 
(for pixels that receive little illumination) is that of a linear increase in the measured flux over time, 
with a slope that varies inversely with the logarithm of the illumination level. 

This effect is important for two reasons. First, it means that the observed 3% increase in 
the measured background flux at 8 /nm over the period of the observations is most likely not the 
result of a real change in the zodiacal background, but is instead another example of this detector 
ramp. Although the increased noise and smaller size of the background at 5.8 ^m obscures this 
effect, there appears to be a similar upward trend. Thus, rather than calculating the background 
in each image individually and subtracting that value, we subtract a constant background of 4.81 
MJy Sr~^ per pixel from all the 8 /im images, and 0.96 MJy Sr~^ per pixel from all 5.8 /im images. 
This background is calculated as the median background value during the last hour of observation, 
when presumably the background is closest to the value it would ultimately obtain. We note that 
this choice has a negligible effect on our final eclipse depths, as the background constitutes only 
0.2% and 1.7% of the signal in our aperture at 5.8 and 8.0 /im, respectively. 

This effect also produces a 1% increase in the measured flux from the star at 8.0 ^m over the 
period of these observations (see Figure 1). There does not appear to be any ramp visible at 5.8 /im; 
instead this data shows a general downward trend. Unlike the detector ramp at low illumination 
levels, the ramp for higher illuminations has an asymptotic shape, with a steeper rise in the first 30 
minutes of observations. We discard the first 30 minutes of data in both the 5.8 and 8.0 /xm channels 
and fit the remaining binned time series from our 3.5 pixel aperture with a quadratic funct ion of 



In At, where At is the change in time from the start of the observations. iKnutson et al.l (l2007al ) used 
this same functional form (with additional degrees of freedom) to describe the detector ramp in the 
8 //m channel over a period of 33 hours of continuous observations of HD 1 89733, and it accuratel y 



captures the behavior of this ramp for a range of illumination levels. Unlike iKnutson et al.l (|2007al ). 
we do not attempt to correct each of the pixels in the images individually for this ramp, as this 
is not necessary for our analysis. Moreover, the lower fluxes, shorter time frame (5 hours instead 
of 33) and reduced cadence of our data (from cycling between the four detectors) make it difficult 
to characterize this effect accurately at the pixel level. Instead we assume the effect of detector 
ramp in the binned flux in our 3.5 pixel aperture will have a shape that is representative of the 
flux-weighted average of the ramp effect for individual pixels. 

We fit both the quadratic function of In At and the theoretical eclipse curve to the data 
simultaneously using a Markov Chain Monte Carlo method as described in ^2.1.11 As before, the 
distribution of values was very close to symmetric in all cases, and there did not appear to be any 
strong correlations between variables. Best-fit eclipse depths and times from these fits are given 
in Table [H and the time series both before and after correcting for detector effects are shown in 
Figures [T] and EJ respectively. 
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2.2. IRS Photometry 

As a check and to ensure a consistent methods between different bandpasses, we re-analyze 
the raw time series gathered with the Spitzer IRS instrument (Houck et al. 2004) at 16 //m (prio r 



to normahzation using the time series for the M dwarf companion) from iDeming et al.l (|2006l ). 
Importantly, these data were obtained on UT 2005 November 17, which falls between the times of 
our observations with MIPS (§2.3) and IRAC (§2.1). We fit the 16 /im time series with the same 
quadratic function of In At as the 5.6 fim and 8.0 /im data. As before, we fit this normalization 
and the eclipse curve to the data simultaneously using a Markov Chain Monte Carlo method with 
the error for individual points set equal to the standard deviation of the points after the end of 
the eclipse, calculated after fitting this section of data with a linear function of time to remove 
the overall slope. This data shows the same type of detector ramp as the 8.0 /im data, and we 
exclude from our fit the first 30 minutes of data where this ramp is steepest. Our estimate of the 
eclipse dept h, 0.519 ifc 0.020%, is slightly shallower and has a smaller uncertainty than the value 



reported bv iDeming et al.l (|2006l ). 0.551 it 0.030%. We also performed a separate analysis of the 
time series assuming a simple quadratic function of time, which we found yielded the same eclipse 
depth and uncertainty as that stated in Deming et al. (2006). We note that the analysis performed 
by those authors differed from ours in that they first divided the time series of HD 189733 by a 
fourth-order polynomial-fit in In At to the time series of the M-dwarf companion, and subsequently 
fit a quadratric function of time to the out-of-eclipse portion of the corrected time series. They 
also reported the results of a linear fit to the corrected time series, for which they obtained an 
eclipse depth of 0.521%, very similar to our final value, stated above. As discussed in ^2.1.21 we feel 
that fitting with a quadratic function of In At represents both a more robust normalization than a 
quadratic function of time, as well as being a better description of the shape of the detector ramp. 
The reduced is 1-07 for the fit that uses a quadratic function of time, and falls to 1.00 when we 
adopt a quadratic function of In At, indicating that the latter function is indeed a better fit to the 
data. Importantly, the difference in eclipse depths between the two normalizations is less than la. 



2.3. MIPS Observations 



We observed HD 189733 for 5 hours spanning a time of secondar y eclipse on UT 20 05 November 
11, and obtained 1548 images using the Spitzer MIPS 24 /im array (jRieke et al.ll2004l ). The MIPS 
observing sequence dithers the target position amongst 14 distinct positions on the detector array 
(see §8.2.1.2.1 of the Spitzer Observer's Manual). Due to the variations in apparent sensitivity at 
each position, we elected to treat the series of images at each position as independent data sets. 
After we produce a photometric time series for each position, we normalize each time series to a 
common scale and subsequently combine them to form a single time series, as described below. 



We estimate the sky background in each image from a 40 x 40 pixel box centered on the star, 
trimming all pixels more than 3a away from the median value and fitting a Gaussian function to 
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the central region of a histogram of the remaining pixels. We find the average background over 
the duration of our observations is 21.9 MJysr"^. We then subtract this background from each 
image individually, and fit the remaining flux in a 10 x 10 pixel subarray (which is large enough 
to encompass the first Airy ring) centered on the approximate position of the star. We found that 
varying the size of the subarray we used did not affect the final time series. We use a MIPS model 
point spread function (PSF) for a 5000 K point soured and allow both the x and y position of the 
PSF as well as a constant scaling factor to vary in our fits. We correct for the effects of transient 
hot pixels in our subarray by making a histogram of fiux values at a given pixel position over the 
set of images (of which there were 108 — 126) at a given nod position and flag 4cr outliers. We then 
assign a statistical weight of zero to these values in the fitting process described above; 95% of our 
images have one or fewer of these hot pixels. We discard images with 5 or more hot pixels in the 
subarray; there are only 10 of these images in the entire time series. We note that the M-dwarf 
companion to HD 189733, is included within our 10 pixel subarray; we give zero statistical weight 
to the values within a 3 x 3 pixel box centered on the position of the companion in our fits. 

We note that the measured background periodically drops by 1.5%, co rresponding to t he firs t 



image in each bundle of images spanning the sequence of nod positions. iDeming et al.l (|2005bl ) 
noted the same feature in their data, and chose to divide the measured fiux by the background in 
order to remove the effect of this drop. This procedure would introduce significant noise in our 
data, because (unlike the case of HD 209458) our target star has a flux twice the value of the 
background. Thus we elected to simply discard the 136 images with these low sky backgrounds. 

For times later than 1.5 hours after the start of observations, the flnal time series at each of 
the 14 nod positions is well-described by a combination of the theoretical eclipse curve and a linear 
function of time. At earlier times, there appears to be a slight upward trend, possibly corresponding 
to the detector ramp described in ^2.1.2[ We fit the entire time series, but find that clipping this 
initial segment of the data in our fits changes the final measured eclipse depth by only 0.5(t. We fit 
a single eclipse curve and 14 linear functions of time (one for each of the 14 nod positions) to the 
data simultaneously using a Markov Chain Monte Carlo method as described in §2.1.11 trimming 
3(T outliers in our fit. A binned version of the combined time series appears in Figures [T] and O and 
we state our estimates of the eclipse depth and its time of center in Table 1. 



3. Discussion and Conclusions 

With the 6 distinct band-integrated fiux ratio measurements in hand, we proceeded to compare 
these values with model predictions. Following Barman et al. (2005), two model atmospheres for 
HD 189733b were constructed, one assuming that the absorbed stellar fiux is confined to the day 
side while the other model uniformly distributes the stellar fiux over the entire planet. Figure [3] 



^Available at |http://ssc. spitzer.caltech.edu/mips/psf.html 
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compares these models and the best-fit Planck curve (corresponding to a temperature of 1292 K) to 
the measured day side secondary eclipse depths. We refer the reader to Barman (2008) for a more 
detailed examination of atmospheric models of HD 189733b in light of our data, but we discuss 
several of the salient points below. 

Our data show clear deviations from the Planck curve. The greatest discrepancy occurs at 
3.6 //m, for which the planet-to-star contrast ratio is more than twice as great as would be expected 
from pure blackbody emission. This emission feature is a prediction that is common to nearly all 
published theoretical models (e.g. Barman et al. 2005; Burrows et al. 2006; Fortney et al. 2006; 
Seager et al. 2005). It arises from a hole in the atmospheric opacities near 4 /xm; significant 
absorption from both water and carbon monoxide in adjacent wavelengths regions causes flux to be 
squeezed out through this wavelength span of relatively low opacity. Furthermore, the data exhibit 
a clear trough from 4.5 — 8.0 /xm, which is consistent with water absorption. This finding is at 
odds with the results from Grillmair et al. (2007), who found that the spectrum of HD 189733b 
was remarkably fiat across the 7.5—14.7 jum bandpass and, in particular, that it did not show 
the downturn at blue wavelengths expected from water absorption. However, the 8.0 and 16 /um 
photometric data in Figure 3 are themselves only modestly inconsistent with the relatively noisy 
spectrum from Grillmair et al. (2007). It may be that while water is indeed present and seen 
in absorption at shorter wavelengths, additional spectral features are present in the 8 — 10 /nm 
region that mask the expected absorption due to water. Fortney & Marley (2007) have noted that 
the Grillmair et al. (2007) result is itself difficult to reconcile with previously published 8.0 /xm 
photometry (Knutson et al. 2007a), since it would require the presence of physically implausible 
opacity sources. 

The presence of water in absorption is in stark contrast to recent findings that water is in 
emission for HD 209458b (Knutson et al. 2008; Burrows et al. 2007a) suggesting fundamental at- 
mospheric differences between these two planets. Several authors (e.g. Fortney ct al. 2008; Burrows 
et al. 2008; see also Hubeny et al. 2003) have discussed how the atmospheres of hot Jupiters might 
bifurcate into two groups, with the presence or absence of a temperature inversion (as deduced 
from the presence or absence of emission features in the spectrum) arising under different levels of 
external irradiation. As noted by these authors, a candidate for the physical origin of this inversion 
is the presence of gaseous TiO and VO in the upper atmospheres of the more strongly irradiated 
planets. Interestingly, HD 209458b has an equilibrium temperature of 1449 it 12 K and hence is 
significantly more strongly illuminated than HD 189733b, which has an equilibrium temperature 
of only 1201 it 12 K (Torres et al. 2008). While it is surely too early to state that the presence of 
a thermal inversion in the atmosphere of HD 209458b and the apparent lack of such an inversion 
in HD 189733b indeed result from their differing levels of irradiation, the connection is intriguing 
and worthy of further study. Moreover, not only do the overall bolometric irradiances of the two 
planets different, but the GOV spectral type of HD 209458 ensures a substantially bluer stellar 
spectral energy distribution than that of the K2V spectral type of HD 189733. We speculate that 
the enhanced short-wavelength irradiation for HD 209458b might generate dramatically different 
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photochemical products and hazes than are present on HD 189733b (see Marley et al. 2007 for a 
discussion of photochemical hazes on hot Jupiters). Since the presence of hazes and their approx- 
imate size distribution may be inferred from transmission spectroscopy (e.g. Pont et al. 2008), we 
urge future modeling efforts to consider the constraints offered by those studies in conjunction with 
the emission spectroscopy results we present here. We note that within the Solar system, planetary 
atmospheres exhibit a wide variety of photochemical products and distinct physical mechanisms 
for temperature inversions, and we should expect no lesser diversity amongst the atmospheres of 
planets orbiting other stars. 

The planets also differ substantially in terms of their physical structure: HD 209458 is larger 
but less massive than HD 189733, such that the densities of the two planets, 0.34 it 0.02 gcm~^ for 
HD 209458b and 0.96 ± 0.08 gcm"''^ for HD 189733 (Torres et al. 2008) differ by nearly a factor of 
3. The larger-than-predicted radii of some, but not all gas-giant exoplanets is one of the greatest 
conundrums in the field today. Many authors have sought to explain the radii of the puffy planets 
by invoking an additional source in the energy budget of those planets (e.g. Bodenheimer et al. 
2003; Showman &; Guillot 2002; but see Burrows et al. 2007b for a dissenting opinion). As the 
atmosphere is the energy gateway for the planet, it is natural to ponder whether there exists a 
physical connection between the observed properties of the atmospheres of these planets and the 
significant differences in their physical structure. The best avenue for progress on the observational 
front is to gather data for a planet that represents an extreme member of one of these two classes. 
We note that TrES-4 (Mandushev et al. 2007) has the lowest density (0.21 ± 0.03 gcm-^) of the 
known transiting exoplanets, and, with an equilibrium temperature of 1785 K it is also one of the 
most strongly irradiated. Hence it may be the ideal target with which to confirm or refute this 
proposed connection. 

The 3.6 — 8.0 fim IRAC observations clearly favor a model that assumes no redistribution of 
extrinsic energy to the night side (upper curve in Figure 3), which would predict very cold night 
side temperatures. This inference is in contradiction to the recent 8 iim light curve measurement 
from Knutson et al. (2007a). Those authors determined that the minimum hemisphere-integrated 
brightness temperature was 973 it 33 K and the maximum hemisphere-integrated brightness tem- 
perature was not much larger, 1212 zt 11 K. This relatively modest day-night contrast indicated 
that the planet was efficiently transporting energy from the day side to the night side. The essence 
of the problem is one of energy budgeting: Our IRAC data indicate a very hot dayside, and the 
Knutson et al. (2007a) values indicate a relatively flat phase curve. If one assumes that the planet 
exhibits a similar lack of a day-night contrast at all wavelengths, the overall emission will exceed 
the energy it absorbs even under the most favorable case (an albedo of 0). Although our 16 and 
24 fim values are in agreement with a model with a small day-night variation (the lower curve in 
Figure 3), the emission at these long wavelengths represents only a small fraction of the bolometric 
flux of the planet and thus cannot, in itself, alleviate this concern. A partial reconciliation may be 
that the degree of energy redistribution is depth-dependent, and hence the emission phase curve of 
the planet is strongly wavelength dependent, such that the planet radiates relatively little energy 
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from its night side at shorter wavelengths. A hint of this effect is provided by a recent upper 
hmit on the day side K-hand fiiix (Barnes et al. 2007) implying much cooler atmospheric tem- 
peratures across the deeper ncar-IR photosphere than predicted by the no-rcdistribution model. 
Alternately, it could be that the planet possesses an internal energy source and hence, like Jupiter, 
simply radiates more energy than it receives. One contender for such an energy source is the tidal 
dissipation of the orbital eccentricity e, which is known to be small but not zero (Knutson et al. 
2007a). However, the predicted dissipation rate (Bodenheimer et al. 2003) based on the Knutson 
et al. (2007a) value is far too small to make this a viable explanation. Barman (2008) considers 
the energy budget of this planet in detail. An additional constraint on the day/night temperature 
contrast may come from the recent detection by Swain et al. (2008a) of methane (which, at colder 
temperatures, becomes the dominant carbon species over CO) in the transmission spectrum of 
HD 189733b. While a cool night side would naturally result in a higher methane abundance, we 
note that non-equilibrium chemistry as recently explored by Cooper & Showman (2006) could also 
lead to a homogenization of the CO and methane abundances across the planet. 

The estimate of the eclipse depth in the Knutson et al. (2007a) 8.0 /xm data was 0.3381 it 
0.0055%, which is 2.3cr less deep that the value we present here (Table 1). This raises the unnerving 
possibility that the dayside emission from the planet may be varying in time, as discussed by 
Rauscher et al. (2007a). This would be both a blessing, in that it would afford us the opportunity to 
study weather on this planet, and a curse in that it would significantly complicate the interpretation 
of broadband spectra assembled from measurements across several secondary eclipses. Star spots 
cannot account for this change: Although the star displays a well-documented periodic brightness 
variation due to star spots (Winn et al. 2007a), we find upon scaling the variation observed at 
visible wavelengths to our 8.0 /um bandpass that star spots would alter the measured eclipse depth 
by only 0.0013%, a negligible effect. 

The observed arrival times of the secondary eclipses provides an exquisite constraint on ecosw, 
where e is the orbital eccentricity and lu is the argument of periastron (see e.g. Charbonneau 2003). 
After accounting for light travel time across the system, Knutson et al. (2007a) found that the 
secondary eclipse of HD 189733 occurred 120 it 24 s later than the prediction for a circular orbit, 
indicating a residual orbital eccentricity of at least 0.10 ± 0.02%. Our data are not of sufficient 
precision to confirm this offset. We find that in 5 of the 6 bandpasses, the eclipse center occurs at a 
time that is consistent with the expectation of a circular orbit (assuming the value of P from Winn 
et al. 2007a and the value of Tc from Knutson et al. 2007a; see Table 1). Intriguingly, the eclipse at 
3.6 fim was observed to arrive 5.6 ± 0.8 minutes later than expected. Since this offset is chromatic, 
is surely isn't due to orbital eccentricity, but instead may arise from strong brightness variations 
across the visible hemisphere of the planet (Williams et al. 2006; Rauscher et al. 2007b). Despite 
the seemingly high significance of this offset, we caution that this band is both the one with the 
largest instrumental corrections arising from intra-pixcl sensitivity variations (Figure 1) and the 
largest residual systematic noise (Figure 2). Nonetheless, it may be that with greater cadence and 
more photons, one could detect the perturbations at ingress and egress that result from a brightness 
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variation across the dayside hemisphere. By electing to observe in a single band and obviating the 
need to continually re-point the telescope as was necessary for the current data, one could realize 
an improvement of more than a factor of 4 in both cadence and total photons. Importantly, such 
observations could be accomplished during the post-cryogenic phase of the Spitzer Mission. 

One of the salient features of the population of hot Jupiter exoplanets is its diversity. These 
objects occupy what might first appear to be a rather limited range of orbital phase space. However, 
much hard work over the past 8 years has found 28 such planets that transit, providing precise 
constraints on their masses and radii. The scatter in their physical properties is astounding, and has 
challenged a community of theorists that must, after all, attempt to unify the bestiary. Similarly, 
with our first detailed comparative studies of the atmospheres of these exoplanets, we are once again 
confronted by the differences amongst objects in seemingly similar environments. This diversity 
is both shocking and invigorating to the theoretical enterprise. We have only one more year to 
gather observations longward of 4.5 fim, at which time Spitzer will exhaust its cryogen and the 
unique window into the atmospheres of these exoplanets will close. Fortunately, we have learned 
that nnich of the key information regarding these atmospheres will be, in principle, contained in 
the two bandpasses afforded by the post-cryogenic phase of the Spitzer Mission. If it is approved, 
the Spitzer Warm Mission will therefore be the pre-eminent tool for investigating the atmospheres 
of these distant worlds until the launch of JWST in the middle of the next decade. 

This work is based on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
Support for this work was provided by NASA through an award issued by JPL/Caltech. HAK was 
supported by a National Science Foundation Graduate Research Fellowship. 

Facilities: Spitzer(IRAC, IRS, MIPS). 
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Fig. 1. — Raw time series observations of HD 189733b spanning predicted times of secondary 
eclipse. The central wavelength of observation for each data set is (from top to bottom) 3.6, 4.5, 
5.8, 8.0, 16, and 24 //m. Each time series is binned in 3.5 minute intervals, normalized, and plotted 
with a distinct constant offset for clarity. The overplotted curves show the best-fit corrections for 
detector effects (see 
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Fig. 2. — Time series observations of HD 189733b after correcting for detector effects (see The 
central wavelength of observation for each data set is (from top to bottom) 3.6, 4.5, 5.8, 8.0, 16, 
and 24 /xm. Each time series is binned in 3.5 minute intervals, normalized, and plotted with a 
distinct constant offset for clarity. The best-fit eclipse curves are overplotted. 
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Fig. 3.— Comparison between IRAC (3.6, 4.5, 5.8, & 8.0 Aim), IRS (16 urn), and MIPS (24 /xm) 

eclipse depths (circles) and model planet-star flux ratios (solid lines). The top black curve is a 
model in which the emission of the absorbed stellar flux is constrained to the day side only. The 
lower black curve is a model with uniform energy redistribution over the entire planet. Integration 
of these synthetic spectra over the Spitzer band-passes are indicated with diamonds. The flux ratio 
expected under the assumption that the planet radiates a Planck spectrum with a temperature of 
1292 K planet is shown as the dashed line. 
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Table 1. Best-Fit Eclipse Depths and Times 



A {fira) 


Relative Eclipse Depth 


Center of Eclipse (HJD) 


O-C (mill.)'' 


3.6 


0.00256 ± 0.00014 


2453699.28340 ± 0.00044 


5.6 ±0.8 


4.5 


0.00214 ± 0.00020 


2453699.27965 ± 0.00087 


0.2 ± 1.3 


5.8 


0.00310 ± 0.00034 


2453699.28281 ± 0.00153 


4.7 ±2.2 


8.0 


0.00391 ± 0.00022 


2453699.28013 ± 0.00102 


0.9 ± 1.5 


16^ 


0.00519 ±0.00020 


2453692.62302 ± 0.00056 


-1.1±0.9 


24 


0.00598 ± 0.00038 


2453685.96845 ± 0.00094 


0.5 ±0.7 



^Our analysis of data from lDeming et al.l ||2006| '). 

^Observed — Calculated times of center of eclipse, with the predicted based on 
the time of center o f transit Tc from lKnutson et al.l ll2007al ) and the period P from 
IWinn et al.] l l2007al '). 



